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SUMMARY 


Analytical expressions are derived for the flapping, the thrust, 
the torque, and. the profile -drag power of a hinged rotor that are 
applicable to high-speed helicopters and to certain types of convertible 
aircraft. The development differs from that used in the standard rotor 
theory in that no limitation is placed on the magnitude of the blade- 
section inflow angles and differs also in the treatment of the reversed- 
velocity region. The equations may be used to calculate the performance 
of a lifting rotor at any angle of attack either directly or, prefer- 
ably, from charts . 


INTRODUCTION 


Present rotor theory (references 1 to 3 ? for example) has proved 
to be entirely adequate for predicting the present-day performance of 
autogiros and helicopters. With the envisioned doubling of the top 
speed of present-day helicopters, however, and with the advent of con- 
vertible aircraft, a review of the assumptions on which the standard 
rotor theory is based was considered desirable in order to determine 
the extent to which the theory could be applied to these improved con- 
figurations. A review of the theory revealed that the premise that 
rotor-blade- sect ion inflow angles 0 are small enough to allow the 
usual small-angle assumptions that cos 0 is equal to unity and sin 0 
is equal to 0 would not apply to the inflow angles generated at the 
rotors of high-speed, high-performance helicopters or to certain types 
of convertible aircraft with rotors which operate through a 90° range 
of angle of attack. Also, for helicopters operating at tip-speed ratios 
close to 0.5, the high in flow angles (and section angles of attack) 
usually associated with high-speed flight, together with the relatively 
large area affected, make the contribution of the reversed-velocity 
region much more significant than it is at the normally low values of 
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tip-speed ratio. Consequently, an investigation was made to remove the 
small-angle and reversed-velocity limitations from the present rotor 
theory. Just as in the present theory, however, no attempt was made in 
this investigation to account for blade stall in that part of the rotor 
disk outside of the reversed-velocity region or for compressibility 
effects on the blade sections. The results of the investigation are 
reported herein. 

A somewhat similar investigation was made at the Georgia Institute 
of Technology under the sponsorship and with the financial assistance 
of the MCA (reference 4) . Both of the investigations made use of the 
idea, first advanced in reference 4, that great simplifications could 
be effected in the equations for rotor characteristics by representing 
the angle of attack of a blade section by its sine. (This idea is 
explained more fully in the section entitled ''Basis of Analysis.") The 
investigations differ, however, in many basic respects, some of major 
and others of minor significance. The most important differences are 
as follows: The analysis reported herein is based on the same reference 

system of axes (that is, the axis of no feathering) as was used in 
previous MCA works on rotating -wing-aircraft theory (see references 1 
and 2 )} whereas the analysis of reference 4 is based on a system with the 
axis perpendicular to the rotor tip-path plane. The present analysis 
also makes use of the "energy" method of calculating rotor losses, a 
method that readily lends itself to the construction of simplified 
performance charts, such as those published in reference 3- The analysis 
of reference 4 uses a "balance of force" method for calculating rotor 
losses that entails the calculation of the longitudinal component of 
force in the tip-path plane (that is, the so-called "H" force) . The two 
investigations also differ significantly in the manner in which the 
reversed-velocity region is handled. It is expected that the analysis 
presented in this paper will be most applicable to aircraft having 
rotors with flapping blades, whereas an analysis based on tip-path-plane 
axes might be applied more conveniently to rigid-rotor aircraft. 


SYMBOLS 

Physical Quantities 
b number of blades per rotor 

R blade radius measured from center of rotation, feet 

r radial distance to blade element, feet 

ratio of blade -element radius to rotor -blade radius (r/R) 


x 
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c blade-section chords feet 


c e equivalent blade chord (on thrust basis), feet 


a rotor solidity (bc e /itR) 

9 blade-section pitch angle; angle between line of zero lift of 

blade section and plane perpendicular to axis of no 
feathering, radians 

9 q blade pitch angle at hub, radians 

difference between hub and tip pitch angles; positive when 
tip angle is larger, radians 

mass moment of inertia of a blade about flapping hinge, slugs 
per square foot 

7 mass constant of rotor blade; expresses ratio of air forces 

to mass forces (cpaR^/lj) 

p mass density of a£r, slugs per cubic foot 


Air -Flow Parameters 

V true airspeed of helicopter along flight path, feet per 

second 

ft rotor angular velocity, radians per second 

a rotor angle of attack; angle between axis of no feathering 

(that is, axis about which there is no cyclic-pitch change) 
and plane perpendicular to flight path, positive when axis 
is pointing rearward, degrees 

v induced inflow velocity at rotor (always positive), feet per 

second 

P tip-speed ratio 
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inflow ratio ^V_gi^_ JL vj 

blade azimuth angle measured from downwind position in 
direction of rotation, radians 

component at blade element of resultant velocity perpendicular 
to blade -span axis and to axis of no feathering, feet per 
second 

nondimens ional component of resultant velocity at blade 
element (Up/£IR) 

component at blade element of resultant velocity perpendicular 
both to blade-span axis and Up, feet per second 

nondimens ional component of resultant velocity at blade 
element (Up/ftRj 

resultant velocity perpendicular to blade -span axis at blade 
element, feet per second 

nondimens ional resultant velocity at blade element (u/ftR) 

inflow angle at blade element in plane perpendicular to blade - 
span axis, radians ^tan -1 

blade-element angle of attack, measured from line of zero lift, 
radians (0+0) 



Aerodynam i c Characteristics 
section lift coefficient , 

average section lift coefficient in reversed-velocity region 
section profile-drag coefficient 


average section profile-drag coefficient in reversed-velocity 
, region 
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5 0 , 8 1> 5 2 


coefficients in power series expressing c^ as a function 


of “r ( c d Q = 5 0 + 8 l a f + 6 2 a ^ 2 ) 


slope of curve of section lift coefficient against section 
angle of attack (radian measure) 


C 1> C 2, C 3,C]J 
and 

Kq , • • • 


coefficients for use in profile -drag torque and power 
expressions 


L 

T 

Q 


lift,, pounds 

rotor thrust, pounds 

rotor-shaft torque, pound-feet 

rotor-shaft profile-drag power, pound-feet per second 


Uip 

C Q 

C *c 

B 


thrust coefficient 


\i® 2 p(nR) 2 y 

rotor-shaft torque coefficient 


Q 


,jft 2 p(flR)%V 

rotor-shaft profile-drag power coefficient 




tip-loss factory hlade elements outboard of radius BR 
are assumed to have profile drag but no lift 


Rotor-Blade Motion 

P . blade flapping angle at particular azimuth position measured 

from plane perpendicular to axis of no feathering, radians 

/ 

Sq constant term in Fourier series that expresses P (radians) 

hence, the rotor coning angle 



coefficient 

of 

cos nty 

in Fourier 

series 

that 

expresses 

P 

8 n 

coefficient 

of 

sin mjr 

in Courier 

series 

that 

expresses 

P 



6 


NACA TN 2665 


Subscripts: 

1 induced 

o , profile 

f "forward-velocity" region 

r reversed-velocity region (except in a^.) 

2 component due to lift vector 

d component due to drag vector 

• 75R at 0-75 radius 

std standard theory 

ext extended theory 


BASIS OF ANALYSIS 


Although the general development of the thrust, torque, and blade- 
flapping equations that is given in the following section is similar to 
that for the standard theory presented in references 1 and 2, the omis- 
sion of the small-angle assumptions in the equations of this paper and 
the consideration of higher tip-speed ratios result in differences 
between the development of the two sets of equations. These departures 
from conventional rotor theory are now discussed. 

(1) The elemental thrust at a blade section used in deriving the 
thrust and thrust-moment equations is taken as the projection of the 
elemental lift vector on the axis of no feathering instead of being 
assumed equal to the elemental lift. The component of the profile-drag 
vector in the thrust direction is also considered in calculating the 
thrust and thrust moment in the reversed-velocity region. (As a result 
of many sample computations, it was concluded that the profile-drag 
component is negligible everywhere except in the reversed-velocity 
region.) Similarly, exact components of the elemental lift and drag 
vectors are used in deriving the torque equations. (See- fig. 1.) 

(2) The resultant velocity U is used instead of lip in computing 

the lift and drag at each element. (See fig. 1.) Radial (that is, 
spanwise) components of velocities are ignored as usual, although such 
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velocity components may significantly affect the blade drag at extremely 
high tip-speed ratios. 

(3) When flight up to tip-speed ratios of 0.5 is considered, the 
contributions of the reversed-velocity region become a significantly 
large part of the total thrust, torque, and power produced by the rotor, 
inasmuch as high inflow angles (and section angles of attack) are 
usually present at the high tip-speed ratios that produce the relatively 
large reversed-velocity areas. (This statement is principally true for 
the helicopter power-on flight condition. For other flight conditions, 
such as partial -power helicopter descents or autorotative operation at 
high tip-speed ratios, large reversed-velocity areas may exist with 
small inflow angles.) The contributidn of the reversed-velocity region 
to the rotor characteristics is computed in the rotor theory of refer- 
ences 1 and 2 by assuming that the blade sections in that region, as 
well as those in the forward-velocity part of the disk, remain unstalled 
no matter how high the angles of attack encountered in that region may 
be. Because the standard rotor theory is also based on the assumption 
that the resultant velocity at a blade section can be represented by 
the tangential component of velocity Up, the assumption of no stall in 
the reversed-velocity region does not normally seriously overestimate 
the thrust and accelerating torque, inasmuch as the overestimation of 
the lift coefficient is somewhat compensated for by the underestimation 
of the dynamic pressure at each blade section. The underestimation of 
the decelerating torque and the profile-drag power is much more signif- 
icant, however, because the effects of the two assumptions are additive 
in those cases. When the normal component of velocity Up is included 
in the resultant velocity at each section, as is done in this paper, and 
when large areas of reversed velocity and high section angles of attack 
are present, the effect of the no-stall assumption results in further 
errors because the thrust and accelerating torque contributions are 
greatly overestimated, and the decelerating torque and profile-drag 
power contributions to -the total thrust, torque, and power are 
underestimated . 

These effects were investigated for a sanqple practical operating 
condition (p = 0.4, C T = 0.00362, 6 q = 24°, = -12°) by numerical 

step-by-step calculations. The results showed that the thrust produced 
by the reversed-velocity region, calculated on the basis of unstalled 
flow and with Up considered in the section resultant velocity, is 
about 23 percent of that produced by the rest of the disk, whereas the 
contribution is only 8 percent when calculated on the basis of stalled 
flow. For the same sample case, the profile-drag power contribution of 
the reversed-velocity region is l8 percent of that produced by the rest 
of the disk when calculated by means of the no- stall assumption, whereas 
this contribution. is 51 percent when calculated with stall considered. 

In addition to illustrating the errors that could be introduced by not 
considering the effects of stall in the reversed-velocity region at high 
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speeds, these data point up the fact that the contributions of the 
reversed-velocity region can be quite significant. 

In order to account for the stalled flow in the reversed-velocity 
region in a relatively simple manner, advantage is taken of the fact 
that the average section angle of attack in this region for most tip- 
speed ratios and power-on flight conditions is extremely and uniformly 
high. (For these conditions, a typical value of section angle of attack 
is 6o°.) The section lift and drag coefficients in this region are 
therefore assumed to be constant. The trust, torque, and power equations 
are then developed by first calculating the contributions of the 
"forward-velocity" region (that is, the area outside the reversed- 
velocity region) and then algebraically adding to them the contributions 
of the reversed-velocity region computed on the basis of constant lift 
and drag coefficients. 

Calculation of rotor characteristics in the reversed-velocity 
region requires certain further assumptions of a semiempirical nature 
even though such assumptions are not necessary for the forward-velocity 
region. The approximations used in the reversed-velocity region are 
pointed out as they arise in the development of the analysis and are 
subsequently shown to be satisfactory over a wide range of operating 
conditions. v 

(4) Preliminary investigations of the problem of computing rotor 
characteristics without making the usual small-angle assumptions for 0 
and 0 showed that the resulting equations for rotor thruBt, torque, 
and blade flapping motion are too unwieldy for analytic solution. In 
order to make the equations amenable to practical solution, advantage 
is taken of the fact that, although the pitch angle or the section 
inflow angle separately may be too large for the small- angle assumptions 
to apply, their sum may not be if the rotor Is to operate with the 
significant elements not highly stalled. It is therefore permissible 
to use Oj, and sin otp interchangeably. (Even at a^. = 20°, the 
error due to this assumption is only about 2 percent.) In practice 
the substitution of sin cl,, for a^. in the expressions for the lift 
and drag coefficient of a blade element makes It possible to develop 
the equations for rotor forces without the previously adopted restric- 
tions of 0 and 0 to small angles. 

The analysis given in the following section is developed, as are 
those given in references 1 and 2, for linearly twisted, rectangular 
blades with the flapping h i nge located on the rotor center and perpen- 
dicular to the rotor axis and 'to the blade span. All velocities, forces, 
and moments are referred to the axis of no feathering. An explanation 
of this system of axes, together with means for applying the equations 
based on this system to pure feathering or combination flapping- 
feathering systems, is given in reference 5- 
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AMIXSIS 


Thrust.- From reference 1 , the velocity components at a rotor 
blade element may be expressed nondimens ionally as 


% 

Ufp = jjg- = x + M. sin \Jr 


( 1 ) 


U P " ?Jf ~ ** + \ ^ a l + ^ a 0 + xb l + \ ^*2) cos t + ^ -xa l + \ M-bg )j sin tJt 


and 


Jj ^ a i + 2xbg| cos 2t + ^ Ubi - 2xa^J sin 2i|r 

U / 2 2^/2 

U = flR = fa + U P ) 


Also, from figure 1, 


or 


or 


Up = U cos 0 

Urp = U COB 0 

Up = U sin 0 

Up = u sin f 


( 2 ) 


( 3 ) 


W 


( 5 ) 


Again, from figure 1, the rotor thrust at a blade element with the 
drag component neglected is found to be 
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dT = dL cos 0 

where, from simple blade -element theory, 

dL = c 7 ipU^c dr 
l 2 


( 6 ) 


(7) 


On the assumption that 


CLj. = sin oLp 


( 8 ) 


and that the lift-curve slope is constant, Cj may be expressed as 


Cj = a sin = a sin^0Q + x0^ + 0) (9) 

If equation (9) is expanded by the standard trigonometric identities, 
and the result combined with equations (k) to (7)> then 

= ipac ["jj3in 0 q cos x0 1 + cos 0 q sin xe^Uj, 2 + 


(cos 0 q cos x0-^ - sin 0 Q sin x0-JUpUpj (10) 

The total thrust produced over the "forward- velocity" region (that 
is, that part of the rotor disk over which the flow Is not reversed) by 
a rotor of b blades is 


T f 


Jb_ 

2n 



dr dt - ~ 
dr 2.n 



it P-pR sin 


dT 

— dr dUr 


( 11 ) 


Before equation (10) is substituted into equation (ll) and the 
indicated integrations are performed, it was found desirable first to 
expand cos x0j_ and sin x0^ by the following terms of the power 
series: 
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COS X 0 -j_ 


1 


M 2 , M h 

21 if! 




sin x0^_ = (x0j^ 



( 12 ) 


The series expansion was used in the integration instead of exactly 
integrating the terms containing cos x0^ and sin x0-j_ in order to 

avoid a discontinuous result in the integrated answer at 0-^ = 0 and 

to avoid the necessity of having to use a very large number of signif- 
icant figures in obtaining a numberical value of the thrust for small 
values of 0p- ' 

Substitut ing equations (l), ( 2 ), ( 10 ), and ( 12 ) into equation (ll) 
and integrating on the basis that the chord c ’is constant gives an 
expression for the thrust produced by the forward-velocity region. The 
expression can be made nondimensional by substituting into it the expres- 
sions for Crp and a that are listed in the section entitled "Symbols." 

When these substitutions are made, the final expression becomes, after 
simplification. 



sin 0 


B- 


+ ^ I - h “ 3 - e i( x f + h “ 3x ) 


♦ - 8 g) ♦ 


0 3 X b 5 + e ^ B? 
e l x 30 + e i 155 


+ COS 0 Q 


\ + \&> Z f + | + n 3 ^ + 


0i + u 2 . © 2^ , bM 3/z 6 , „2 

0 ilT + 11 T “ j + 0 i [m " x T ) ~ e i [3S + * TB) 


(13) 


The simplifications expressed in equation ( 13 ) were made by neglecting 
all terms in q. of higher order than the fourth and by omitting all 
terms that contributed but negligibly to the final result} the latter 
terms being deter min ed from numerical examples representing extreme 
flight conditions. These same simplifications were also made in the 
final expression for the flapping, rotor torque, and rotor profile- 
drag power coefficients which are subsequently developed. 

The thrust produced by the lift vectors in the reversed-velocity 
region T r ^ is calculated on the basis that the lift coefficient remains 
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c onst ant , so that 



, 1 _ i 

- b ^ cc i s; 



'-pRsinty 


dijr 


U^cos 0 dr 


(14) 


Inasmuch as exact evaluation of the last integral in equation (14) is 
not feasible, an approximation of the U 2 cos 0 term by an expression 
that is readily integrated is necessary. The semi empirical relation- 
ships employed are as follows: 

U 2 COS 0 = UrpU = Upfup - 





(15) 


In the preceeding equation. Up is replaced by the expression 


Up = UpflR = 



(16) 


which is based on the assumption that Up varies linearly from a value 

of X at the blade root to zero at the blade tip. The same expression 
for Up is also used in deriving the reversed-velocity contributions 

to the rotor torque and power. The bracketed expression in equation (15) 
is made a positive quantity by placing absolute bars about X, inasmuch 
as the expression replaces the resultant velocity U, which is always 
positive. 

Evaluating equation (l4) with the aid of equation (15) and making 
the answer nondimens ional gives the expression for the reversed-velocity 
thrust as 



1 2 

13 “ 


M(i - |) 


9* 


(IT) 


It should be noted that c ^ is positive when X is positive, and 
negative when X is negative, inasmuch as the reversed-velocity 
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thrust acts upward in normal power-off flight and downward in power-on 
flight . 

As mentioned in the previous section, the contribution of the 
profile drag to the total thrust is negligible in the forward-velocity 
region. In the region of reversed flow, however, the profile-drag 
contribution T r ^ becomes significant because of large inflow angles 

and high profile-drag values resulting from the stalled flow and is 
given by the following expression: 


1 if 2 * P-HRsint 

T r d = b 2 pcc d a + / train 0 dr (l8) 

0 Jit JO 


The use of the same assumptions employed for equation (15) yields 



Substituting equation (19) into (l8), integrating, and making the result 
nondimens ional gives the following expression 



( 20 ) 


The total rotor thrust coefficient can then be obtained as the sum 
of equations (13), (17)> and (20); thus. 


2C t 
a a 




( 21 ) 


Blade motion .- The motion of a flapping blade can be described by 
a Fourier series that expresses the relation between the flapping 
angle £S and the blade azimuth angle t . To a degree of accuracy 
sufficient for most purposes, the series can be written as 
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3 = Sq - a^ cob t sin \{r - cos 2\|r - bgsin 2 (22) 

Equations for determining the five unknown coefficients in equation (22) 
are obtained by equating the sum of the blade thrust, centrifugal force, 
and inertia moments about the blade -flapping hinge to zero. 

The thrust moment at any value of i|r may be written as 



] 2 Jt 

is used to indicate that the expression enters 

into the thrust moment only in the interval between i]r = n and ^ = 2jt. 
The elemental thrust moment dTr can be obtained by simply multiplying 
equation (10) by r. The last integral in equation (23), which repre- 
sents the moments contributed by the lift and profile-drag components 
in the reversed- velocity region, can be evaluated by means of equa- 
tions (15) and (19)> (As in equation (17), the sign of c^ is dependent 
upon the sign of X.) 


Evaluating equation (23) and following the method used in refer- 
ence 1 results in the following five simultaneous equations for use in 
determining the five flapping coefficients: 



o 
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0 = cos 6q 


li ^ 1 v ^4 




“3~ ^ + £ ^l® + — g - 0.050a2(i 


,4 *2 9i 2 / B 5 


6 2 -5 


-^T a 0 + 


b x B 


sin 0 


0 


e i(- “0 ♦ m % + >\ |§ 


MB ap | o ( 

— +e i 3 -^, 


(25) 


0 = sin 0 


f 4B3 + 0.026** - 8l U y - *1 %■ + g a lE 3 - £ „,/ - 


5 .2 


0.0265^) - ^IB 5 ) - 0 3 fe ? 7 _ ^„ a l B ^ 

j 2\5 ) 1 \ 42 120 / 


+ cos Or 


H ^B 2 
2 


H 2 a-,B 2 |iB3b p , . , o / ^lA a 2/ ], 

a i T + 8 A - i - A - h (* - 4-U £ 


6 32 


16 


2 / 2 k 


4 + ^ !s) - 4 { 4 ) * ^(- =#) 


1 l ^7 

0.0132^ — - 


0.2123A|l|n 2 (l - |) 2 ^ 


£ -bA'4'^- 


32 X ( 1 ■ 2) 1 


3 Cd o 


(2 6 ) 
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*2 ~ F i 


sin 0 


0 


iS 2 + i _ 


6h 


H - 0 n |ia. 


B 




5\ © 


1 (1 


ipi ^ 2 T/ + — la 


b\ 


0 


3“ a l + — 2 ^ + 0-0265XH 3 + O.OlSOa-jH 1 *- + © 1 B^j - 

H T a i + T b 2) - SS ^ - gf E a 0 + 


0-0265 3 

a H 


w( l -ih + ^ + -iK 


(27) 


b = 4 


cos 0 


2 B 2 p.B^ 

T + 3- b i 


B 1 1 

2 “I 16 


0 


1 / M-B" 


2 \ 5 




- sin 0 


0 


5 ,,2 

3 l( a 0 ~6~ + 


“ h i T - f ^ 


MaB) 


Accelerating torque.- As was done in reference 2, the torque arising 
from the inclination of the lift vectors relative to the plane perpen- 
dicular to the axis of no feathering is called the "accelerating" torque 
in the present paper, even though such torque tends to accelerate the 
rotor only in the case of power-off operation. 

The differential expression for the accelerating torque (see fig. 1) 
may he written as 


12 , 

= TjPU c 1 sin 0 cr dr 


( 29 ) 
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which, after substitution from equations (4), (5)> and (9) , becomes 

~ = ^pac|sin^e 0 + ©px) U p U T + cob (8 q + e^Up^jr (30) 


dQj 

dr 


The accelerating torque produced by the forward velocity region is 


then 


Ql f 2* 


dt 


HBR dQ ± b 

_ — (jp - — — 

dr 2* 


P2jt p-|i Rsimlf dQ i 


dt 


Jo 


dr 


dr (31) 


Evaluating equation (31) with the aid of equations (l), (2), (12), 
and (30), and nondimensionalizing and simplifying the result gives the 
following expression 



where 

Cp = X. 2 + |iapX + | H 2 a Q 2 + ^ ^ 2 a x 2 + | ^ 2 bp 2 - | (33) 


( 34 ) 
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c 3 - I ^ + I * <*£ 


(35) 


The accelerating torque produced "by the reversed-velocity region is 



r M-RsimJf 


ipcU^sin 0c ^r dr 


(36) 


The approximation is made that sin 0 = Up^. Using the value of Up 

as given by equation ( 16 ), substituting into equation ( 36 ), integrating 
and nondimensionalizing the answer yields the following expression for 
the accelerating torque in the reversed-velocity region 


( 2 C 


Qi' 


0 a 


= - S l 1 - 1) ! 


(37) 


The total accelerating torque is then the sum of equations (32) 
and (37) 


2C 


% psi 


2C 


Qil 


ca 


0 a 


0a 


(38) 


Decelerating torque .- The decelerating or profile-drag torque pro- 
duced by b blade elements (see fig. l) is 


dQ Q = bc d ipcl^cos 0 r dr 
u u o ^ 


(39) 


where, as in reference 2, c d is expressed as 


% = 5 0 + 6 l a r + 5 2 a r 2 

/ 


m 
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When equation (39) is expanded with the aid of equations (4), (5), (9)> 
and (40), the following expression is obtained: 


d ^O = | 5 0 ' co s + S l^ Bin ( e 0 + + COS^e 0 + 0-pXjUpUfj, 


8 2 |sin 2 ^ 0 + + cos 2 (0 Q + 0-jX^ 

U p U T sin 2(0 q + 0-Lxjj cos 0 


> r dr 


(4l) 


Although the assumption that = sin ctj, was successful in 

eliminating terms involving sin 0 and cos 0 from the integral 
equations for thrust, thrust moment, and accelerating torque, cos 0 
terms still remain in equation (4l). As will he shown in the succeeding 
section entitled "Discussion," a subsequent investigation showed that 
cos 0 could he assumed equal to unity in equation (4l) with small error 
in the final answer for Q 0 in most cases, even in the most extreme 
operating conditions when 0 is very large. The validity of such an 
assumption is aided hy the fact that, in the normal operating angle- of - 
attack range, the 5 q and S 2 terms are of the same order of magnitude. 
Since cos 0 appears in the denominator of the 8 q term and in the 
numerator of the 8g term, the effects of cos 0 tend to cancel each 
other . 


The decelerating torque contributed hy the forward velocity region, 
as represented hy the nondime ns ional coefficient, is 



dijr 


1.0 dCQ 
dx 




''-psinijr dC 


Qq 


dx 


dx 


(42) 


Substitut ing equations (l), (2), ( 12), and the nondimens ional form of 
equation (4l) into equation (42) and integrating results in the 



following expression for 


f 
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where Cj, C 2 j and C 3 are as defined in formulas (33) to (35) and 


0 

•f=- 

II 

l n 2 \ 2 + | yi\\ + ^ »\ 2 

m 

K l = 

8 q + Bi sin d Q + 5 2 sin 2 0 o 

(45) 

V 

5 1 cos 9q + 5 2 sin 29q 


k 8 = 

^-8^ sin 9q + 2&2 cos 20 q)0-j_ 

(47) 
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% = ^ 5 i cos e o ~ 25 2 sin 20 O^ e l 2 (48) 

K 10 = 5 1 sin e 0 - | S 2 cos 2e o) e i 3 (4-9) 

K 11 = 5 1 cos S 0 + f S 2 sin 20 O ) 9 l k (50) 

K2 = Kf ! ' (51) 

K3 = | ( 52 ) 

K i^ = | % (53) 

% = % K 10 e i (54) 

k 6 = 5 K u e i (55) 

Kj2 = ^20 ®0 + 15 cos 2 ®o)®l^ ( 56 ) 

k 13 = 5 2 cos 2 0 o ( 57 ) 

K l4 = _5 2 sin 20 O 6 1 (58) 


Ki 5 = -8 2 cos 2e 0 e 1 2 


(59) 
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Ki 6 = | & 2 sin 20Q0-L3 (60) 

K 17 = i 5 2 cos 20 q 0^ (6l) 


K^g = - A &2 ■ sin 20Q0-L 5 (62) 


It should be noted that the , K constants defined by the preceding 
equations are functions only of the pitch angle, blade twist, and the 
section drag characteristics. 

The decelerating torque coefficient for the reversed-velocity 
region is 


= k 



pRsin'lr 


Jpbcc d U 2 cos 0 r dr 


(63) 


If the value of U 2 cos 0 as given by equation (15) is substituted into 
equation (63), made nondimens ional, and then integrated, the result is 

0 ), ■ - 1 ) * A 

The resultant decelerating torque contribution is obtained as the 
sum of equations (43) and (64): 


(64) 



(65) 


Profile-drag power .- The differential expression for the profile- 
drag power can be written as 


dP Q = 


1 3 

b gPcc d U dr 


( 66 ) 
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Expanding equation ( 66 ) with the aid of equations (3), (4) ( 5 ) ( 9 ) 

and (40), and expressing the results in nondimens ional form, gives ’ 



0 

2 






Up 

cos 0 



+ 0 x xj 


+ 


^Up 3 + Up^urpj sin^ 0 + e^xj 


+ 6 r 


^T 2 p 
U P cos 


( 9 0 + 6 1 X ) 


U T 3 tu 2 

? sln To + e i*) + ip sin s(e 0 + a A 


COS 


(67) 


o As in equation ( 41 ), it was found that the omission of the cos 0 terms 
in equation (67) results in an error of but a few percent in the final 


The total profile-drag power contributed by the foward- velocity 
region can be written as J 





'-psimjr dCp 


dx 


dx 


( 68 ) 


After substitution, integration, and simplification, the final expres 
sion for Cp ^ becomes: 



2k 
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(i c Po) f ■ + ^ 2 ) + k 4 + H + k3 & + H + 4 + ^ 2 ) + 

+ %(| + a2 ) + Jr v 4 + ( x ‘ I ^(i "t + 


s v 


1 " • 4 £, + i K,„ + i K nl + g Kjgj + |§ A + | ^ 3a lj( K 7 + 


t K 8 - 5 "9 + 6 K 10 + 7 Ku - 


I % + 5 % + 5 K 10 + 5 K n + ? K ia) - I; ^7 - s ^ + 

t * 2 ^ 1 - 5 l* 3 <> 0 > 1 ) ( K 19 + | K 20 + | K 21 + f K 22 + 5 K 23 + 

I K 2lt) + + 2 + 5 '"V + 8 112,312 + 2 ^Vs) (2 Kl 5 + 

f Kao + f %L + 5 %> + 5 *23 + T K 2b) * (- ,la ° bl + i ,1! ‘ lb2 ‘ 

I ^l a 2 )(| K 19 + 4 *20 + 5 %. + 5 ^2 + 7 *23 . + 8 **) + 

(l a l 2 + | b l £ + ea 2 2 + a 2 2 ) (5 K 19 + 5 *20 + I K£1 + 7 ^ + 


g % + 5 **) - (l K is)(5 ^ + 5 “V + K “V’ 



(69) 
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where Kp to Kpg are expressed by equations (45) to (62) and 

2 

K-j^ = Sq + 5p sin + 5g cos 9 q 
KgQ = (&± cos 0 q - S 2 sin 20 q^Sj_ 

K 21 = g 5 1 Sln 6 0 ” 5 2 cos 29 0 j 6 l 2 
Kg 2 - f- g &l cos Sq + f 5 2 sin 20 O^ 0 1 3 

K 23 = 0% &1 BjM d 0 + J 5 2 cos 

*24 = (llo 5 1 COS S 0 ' B 82 8±n 20 o) 0 l 5 


(70) 

(71) 

(72) 

(73) 

(74) 
> 

(75) 


The profile-drag power for the reversed-velocity region is 


, • P2« p-pRsin|r 

( F °)r~~^l "J> ***< 


- U 3 dr 


(7 6) 


The variable part of the integral in equation (j 6 ) can be handled as 
follows : 


U- 


= (ut 2 + Up 2 )u = (up 2 + Up 2 ) Up - ||x|ar(i - |jj (77) 


Substituting the value for Up obtained from equation (16) into 
equation ( 77) ^ substituting the result into equation (76 ) , integrating, 
and nondimensionalizing permits the power loss in the reversed-velocity 
region to be expressed by 



( 78 ) 



t 
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The total profile-drag power coefficient is then determined as the 
sum of equations (69) and (78): 



(79) 


DISCUSSION 


Check of simplifications and assumptions .- A number of assumptions 
and simplifications (among which are the assumption that section lift 
and drag coefficients in the reversed-velocity region could be repre- 
sented by uniform stalled values, the simplifications used in the 
reversed-velocity region, and the neglect of some cos 0 terms in the 
profile-drag torque and power equations) have been made in deriving the 
rotor equations. 

In order to determine the effect of these approximations on the 
accuracy of the theory derived in this paper, rotor characteristics 
for a number of different flight conditions were calculated by the 
theoretical equations, and the answers have been compared with those 
obtained by a step-by-step integration performed on automatic computing 
machines. The step-by-step method eliminated all small-angle assump- 
tions and the assumption that the resultant velocity at the blade 
element U (without the radial component of flow) could be replaced by 
the tangential component Urp. The method also permitted section Btall 
in the reversed-velocity region to be considered, although stall was 
not considered by either method in the forward-velocity region. The 
same lift and drag curves were used by both methods in the forward- 
velocity region and were represented by 


•\ 

c- L = 5-73oh 

> 


c d = 0.0087 - 0.02160^. + O.lj-OOcXr 2 


(80) 


In the reversed-velocity region, the analytical method assumed that 


c l 


1.2 



= 1.1 


(81) 
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whereas the step-by-step method utilized the lift and drag Curves shown 
in figure 2. The curves in figure 2 are arbitrary and were devised to 
simulate, in an approximate manner, the characteristics of an airfoil 
section up to high angles of attack. 

The ratios of values of thrust, accelerating torque, decelerating 
torque, and power obtained by the two methods are plotted against the 
inflow angle 0 at the tip of the advancing blade (fig. 3)* The values 
of rotor angle of attack, tip-speed and inflow ratios, and blade pitch 
and twist angles corresponding to each of the flight conditions are also 
indicated in the figure. It can be seen that good agreement exists 
between the extended theory and step-by-step methods of computing rotor 
characteristics throughout the range of inflow angles investigated. 

(The highest values of inflow angles covered correspond to operation at 
tip-speed ratios of 0.5 or to operation at large rotor angles of attack 
such as might occur at steep climb angles. The theory is expected to be 
applicable for any value of 0 encountered in helicopter or convertible- 
aircraft flight.) It may therefore be concluded that the assumptions 
and approximations used in deriving the present theory are justified by 
the simplifications of the relationships they permit and by the good 
agreement of the results of the theory with those obtained by an exact 
step-by-step analysis. 

Comparison between standard and extended theories . - Essentially, 
the difference between the standard and extended theories lies in the 
fact that the extended theory omits the small-angle assumptions in the 
forward- and reversed-velocity regions and takes stall and the profile- 
drag contribution to the thrust into account in the reversed-velocity 
area. For relatively low-speed and cruising level flight, or at moderate 
rates of climb and descent, the two theories should give substantially 
the same results, whereas in high-speed flight or at large angles of 
climb or descent different results would be expected. It Is of interest 
to determine the range of flight conditions at which the two theories 
begin to diverge, and also the maximum extent to which the differences 
in small-angle and reversed-velocity assumptions affect the rotor 
equations . 

A comparison of the thrust equation as herein developed with that 
given by the standard theory in reference 1 reveals that, aside from 
differences In the reversed-velocity terms, the removal of the small- 
angle assumption resulted in the replacing of 9 by sin 9 in the 
equation and in modifying all other terms by cos 9. For a constant 
forward speed (or constant M-) condition, therefore, the blade pitch 
angle 9 lends itself for use as a convenient parameter by which the 
two theories can be compared. 

In order to illustrate the comparison with a sample case, calcula- 
tions by both theories of thrust, accelerating and decelerating torques. 
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and profile -drag power of a sample rotor having "blades with -8° twist 
are plotted in figure 4. The results, which are presented in each 
figure in the form of ratios of values as computed by the two theories, 
are plotted against the three-quarter -radius pitch angle. (The accel- 
erating torque curves of figure 4 are plotted in the form of absolute 
values, rather than ratios, because in the region of low or zero torque 
the ratios would become meaningless.) The maximum values of 0 

shown on the plots would be about the maximum values that might be 
expected in normal helicopter operation. Very much larger angles might 
possibly be employed in the rotors of certain types of convertible air- 
craft with rotors which are also used as propellers. In such cases, 
low values of p may be combined with very high values of 0. Curves 
are shown for three values of tip-speed ratio (0.1, 0-3> and 0-5) and 
for two values of 2Cij/aa (0.03 and 0.01). The higher value of 2Crp/cra 

represents an upper limit from the standpoint of rotor stall at the 
higher tip-speed ratios, whereas the lower value represents the value 
near which most high-speed helicopters will tend to operate. 

A study of the plots of figure 4 indicates that the two theories 
give essentially the same answers (within 10 percent) at low values 
of p over the range of 0.75R values shown. If curves for the low 

values of p were extrapolated to higher 0 values, significant 

differences could be seen to occur at those pitch values at which the 
rotors of certain types of convertible aircraft might be expected to 
operate. At high values of p, the two theories differ by very large 
amounts, the differences varying markedly with 0. Thus the extended 
theory would seem to be more applicable than the standard theory for 
calculating the performance of rotors operating in the power-on condi- 
tion at very large angles of attack or at tip-speed ratios exceeding 
approximately 0.25. 

It might be pointed out that the thrust equation (equation (13)) 
is particularly sensitive to the effects of 0 in that small errors 
arising from the replacement of sin 0 by 0 and cos 0 by 1 are 
magnified in the equation because the sin 0 and cos 0 terms are of 
approximately equal magnitudes. For the power-on case, the two sets of 
terms are subtracted from one another, and the answer, being a small 
difference between two large numbers, can be very much in error. 

Selection of average lift and drag coefficients in re versed- velocity 
region . - As previously pointed out, it is possible, in autorotation for 
example, for large reversed-velocity regions to exist in which the mean 
angle of attack is quite small as compared with the high, usually 
stalled, values in the power-on condition. The proper choice of the 
average lift and drag coefficient values to be used in calculating the 
reversed-velocity region would therefore depend on the flight condition. 
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For most power-on flight conditions, the numerical calculations dis- 
cussed in this section indicate the following values to he satisfactory: 
c 7 = 1-2} c, =1.1. For the autorotative flight condition, values 

of c z and in the neighborhood of 0.5 and 0.1, respectively, seem 

to be reasonable. For other flight conditions, other values of c z and 
c d might be more appropriate. 


LIMITATIONS OF THEORY 


As stated in the preceding section, the extended theory of this 
paper is expected to be applicable for any range of inflow an g les 
encountered in helicopter or convertible -aircraft flight. Also pointed 
out was the fact that the semiempirical approximations used in handling 
the reversed-velocity region were satisfactory over the wide range of 
flight conditions that were investigated numerically, which included 
tip-speed ratios up to 0.5- Although the numerical check of the theory 
tends to inspire confidence in its application to somewhat greater values 
of p, the use of the reversed-velocity region approximations, as exem- 
plified by equations (15)., (l6), and (77), and the omission of the radial- 
velocity component, probably will result in a loss in accuracy of the 
theory at the higher tip-speed ratios, particularly in the calculation of 
the rotor profile -drag torque and profile -drag power. It is difficult, 
in fact, to evaluate accurately by analytical means the degree of error 
in estimating the rotor drag brought about by ignori n g the radial -velocity 
components of the resultant velocities at each of the blade elements at 
any tip-speed ratio. 

The possible reduction in accuracy at tip-speed ratios much beyond 
0.5 does not place a significant limitation on the theory, inasmuch as 
rotors operating at the very high values of p will necessarily be 
lightly loaded, at which condition the relative importance of the profile- 
drag terms becomes less as compared with the other power absorbing items, 
such as parasite and possibly auxiliary -wing drag. If a relatively high 
degree of accuracy is desired at the higher values, the errors introduced 
into the theory by the approximations and omissions should be analytically 

or experimentally investigated. 

Inasmuch as blade stall was considered specifically only in the 
reversed-velocity region, the extended theory, just as standard rotor 
theory, is limited to conditions of operation wherein excessive stall 
does not exist outside the reversed-velocity boundary. The effects of 
moderate amounts of stall on profile-drag power may be estimated by the 
method suggested in reference 6. Blade stall limitations, which have 
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■been shown to occur simultaneously with conditions of optimum perform- 
ance in reference 6, may "be estimated "by the methods given in refer- 
ences 2 and 3 for the standard rotor theory. 

Limited experimental tests and calculations indicate that compres- 
sibility effects on rotors operating at high tip speeds may substantially 
alter the aerodynamic characteristics of the rotor, the primary effect 
being an increase in profile -drag power. The theory of this paper does 
not account for these losses because of the added complexity involved and 
the lack of suitable airfoil data at high Mach numbers and high angles 
of attack. When such section data, as well as more extensive experi- 
mental rotor data, do become available, it is expected that compres- 
sibility power losses may be accounted for in a manner similar to that 
used for stall effects. 

Since some high-speed helicopters and convertible aircraft will 
probably be operating at relatively high disk loading and torque values, 
an investigation of the effect on the equations of the approximation of 
tip losses by a constant tip-loss factor B and of the ommission of a 
rotational -inflow factor was considered advisable. Therefore a compari- 
son was made of the results of the extended theory, modified for the 
vertical -flight condition, with those of the Goldstein propeller theory, 
for a few different vertical-flight conditions. The cases investigated 
included conditions of relatively high inflow, torque, and disk loading, 
as might be encountered by rotors used in certain types of convertible 
aircraft, as well as more conventional helicopter loadings. The com- 
parison indicated that, at least for the vertical-flight condition, the 
use of the widely used constant tip-loss factor B = 0.97 yielded 
answers for rotor characteristics that are in close agreement with those 
given by the more exact Goldstein propeller theory. A similar check for 
the forward flight condition could not be made because of the lack of a 
rigorous method of calculating tip losses in that condition. The ques- 
tion becomes much less important in high-speed flight, however, inasmuch 
as induced losses are comparatively low at high speeds . 


APPLICATION OF THEORY 


The equations of blade flapping, thrust, torque, and profile -drag • 
power were developed in terms of the three fundamental variables; 
namely, X, 0, and p. Different combinations of these variables may 
be ins erted into the equations in order to obtain the characteristics 
of rotors in different flight conditions. Inasmuch as the equations are 
lengthy when applied to twisted blades, if a large number of different 
rotors or flight conditions are to be investigated, the equations can be 
used to prepare charts (similar to those given in reference 3) from which 
the rotor characteristics can be obtained directly. 

Although the equations of this paper were developed for rectangular 
blades, experience (obtained with the theory of reference 2, which was 
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also developed for constant -chord blades) has shown that such equations 
may be applied successfully to blades of conventional tapered plan form 
by the use of an equivalent chord c e for calculating the solidity ratio. 


CONCLUDING REMARKS 


An extension of rotor theory has been presented to obtain analytical 
expressions for rotor flapping coefficients, thrust, torque, and profile- 
drag power that are not limited by the small-angle assumptions for blade 
pitch and inflow angles used in standard rotor theory. In addition, 
because the contribution of the reversed-velocity region becomes signif- 
icantly large in high-speed rotor flight, it was necessary in the deriva- 
tion to account for blade stall in this region. The equations are 
therefore expected to be applicable to high-speed helicopters and to 
certain types of convertible aircraft. Just as in standard rotor theory, 
the equations are limited to conditions of moderate amounts of stall on 
the blade sections outside of the reversed-velocity region and do not 
include Mach number effects. 

A comparison of the equations of the extended theory with those of 
the standard rotor theory indicates that the extended theory should be 
used in predicting the characteristics of rotors operating in the power- 
on condition at very high angles of attack or at tip-speed ratios 
exceeding approximately 0.25- 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya., December h, 1951 
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Figure 2.- Lift and drag curves used in numerical evaluation of reversed- 

-velocity region. 
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Figure 3.- Ratio of rotor character i a tics as computed "by extended rotor 
theory of present paper and an exact numerical integration analysis 
plotted against Inflow angle at tip of advancing blade (l.e., at 
x = 1.0 and + =* 90°) • 
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(a) p = 0.10. (b) \l o 0 . 30 . (c) p = 0.50. 

Figure 4.- Comparison between standard and extended theory over a repre- 
sentative range of flight conditions} 0 = -8°. 
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